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Variational quantum Monte Carlo calculations are reported for the bulk GaAs semiconductor 
in order to present values for the ground-state energy, the lattice constant, the bulk modulus, 
and some derived properties. The statistical accuracy is significantly higher than the remaining 
differences to the experimental values, especially for the total-energy upper bound. The agreement 
with experiment is satisfactory. The results are also compared with those of density functional 
calculations. The accuracy of our results is comparable to the best of those calculations. 

ON \ PACS Numbers: 71.15Nc, 61.50Lt 
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I. INTRODUCTION 



In comparison with the widely accepted and extensively applied local density functional technique (LDA), the 
quantum. Monte Carlo methods (QMC) represent rather new developments for the ground-state determination of 
solids.Eru Especially the variational quantum Monte Carlo method (VQMC) is attractive, because the quality of 
the approximation is controlled by the property of yielding an upper bound, i.e., by the simple rule: the lower the 
energy the better the wave function. The statistical error of an expectation value can be safely estimated. It behaves 
favorably near an eigenvalue where the variance should vanish. Thus this method is really ab initio, not relying on 
any uncontrolled quantity, such as an unknown correlation functional in the case of LDA. This aspect seems to be 
actually important, because tha-detelled results of the density functional theory (DFT) in fact depend sensitively 
on the choice of this functional. E2Tt3 Furthermore, in the past the role of Monte Carlo simulations for homogeneous 
0^ ' systems was to determine the correlation energy to be adapted as functional for inhomogeneous LDA calculations.!! 
"^h , Thus today's QMC calculations on inhomogeneous systems may also be valued by the DFT as it may lead to an 
improved functional. n 

The first paper to apply VQMC to a solid of nonlight elements was presented by Fahy, Wang, and LouieO dealing 
with the ground-state determination of diamond, graphite, and silicon. They proved that such a calculation is feasible 
if one relics, on a pseudopotential to replace the inner shells' contribution. Different from the low atomic number 
systemsjjolj an uncertainty thereby enters into this technique by the construction of the pseudopotential itself. In 
order to trea-Uthe system on an equal footing with LDA, it is adequate to use the nonlocal ab initio potential developed 
by HamannJlil 

In this paper we apply the VQMC method to the comnound semiconductor GaAs to calculate the ground-state 
properties. We closely rely on our former paper on lithiumQ in treating the solid by simulating a finite system and in 
treating various other calculational details. Atomic simulations are carried through to obtain the binding energy with 
respect to these reference values. Hartree-Fock calculations have been peformed to extract the correlation energy. 
The inner shell electrons are included within a nonlocal pseudopotential for the valence electrons. 

In Sec. we briefly expose the theoretical method used by VQMC. The results are discussed in Sec. III. 



II. METHOD 



The Monte Carlo algorithm calculates expectation values of the form 

= J f(x)p(x)dx . 
j p(x)dx ' 

where the integration variable x denotes typically several hundred single variables. The error is of purely statistical 
nature and depends on the number of simulation points as well as on the variance and the autocorrelation function 
of /. The quantum mechanical expectation value of an observable B is calculated by setting 
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(2) 



The smaller the variance of (Bif;)/ip, the smaller the statistical error of the calculation, which can be taken as a 
measure for the many-body wave function ip to be an eigenfunction of B. 

In the present approach we simulated N valence electrons in a supercell containing a limited number of crystal unit 
cells. We sampled systems of 4 and 32 GaAs unit cells to test for convergence. Only the results for the large system 
with 256 electrons are described here. To model the solid, periodic boundary conditions are assumed, which affect 
the Hamiltonian and the wave function as well. The Hamiltonian 

H = T + E c + E\ oc + E n \ oc (3) 

consists of the kinetic energy T, the Coulomb energy E e , and the local part E\ oc as well as the nonlocal part E n \ oc of the. 
additional electron-ion pseudopotential, which is a generalized norm-conserving ab initio atomic pseudopotential.O 
The Coulomb energy includes the nucleon-nucleon, the electron-ion, and the electron-electron part. It is calculated 
by means of the Ewald summation technique.Ej The local part of the pseudopotential can be decomposed into a sum 
over atoms a at site R Q , 

E loc = H>c,a(r?) (4) 

a,R. Q ,s,i 

with rf being the position relative to R a of the ith electron of spin s. The respective sum for the nonlocal pseudopo- 
tential contains a summation over angular momenta I, upon which the wave function is projected: 



x ^ v , , Y lo (0,0)f r , =r Y l ttn r ,M...,ri = r', . . .)dfi r 
L> VaAT) V(---,rf =r,...) 



VWca (r? = r) = V a , t (r) ^ ' Jr =r .,7 Js _ Z \ ■ (5) 
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The spherical angle Q r ' is referred to the direction of r in an arbitrary but fixed polar coordinate system. The 
projection operator demanding high numerical efforts needs only to be calculated where the radial pseudopotential 
V at i(rl) is nonzero, i.e., when r| lies within the radius of the atomic pseudopotential. The integration is done as a 
sum over a grid of four points on the sphere that is exact for I < 2.E£l 

The wave function consists of a product of two Slater determinants and a Jastrow factor: 

<&(]-!, ...,r N ) = D T (r l5 . . . ,r A r /2 )L> i (r Ar/2+ i, . . . ,rjv)exp I -^u^) I . (6) 
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The exponent of the Jastrow factor is chosen in the usual form for solids: 



U (r) = 4(l-e- r/ir ), (7) 



r 



where A is a variational parameter. For a given A the value of F is fixed by the cusp condition, which removes 
the Coulomb singularity. The Slater determinants D^,D^ consist of parametrized hybrid bonds in the tetrahedral 
directions f from the arsenic ion 

0(r) = < y s b +/3^ y a b , (8) 

< S b(r) = 7 A >f[(r - Ra s )/C As ] + E ^< S t( r - R As)/C As L (9) 

i 

CbM = 7 Ga ^? a [(r - R Ga )/C Ga ] - £ n^[(r - R Ga )/C p Ga ], (10) 

i 

where f3, 7, and C a- re variational parameters and i runs over the Cartesian components [x, y, z). The quantities £ 
describe the contraction of the wave functions as usual for atomic orbitals in solids and reduce the additional extent 
introduced by the Jastrow factor. The coefficients 7 and (3 represent the adjustable parameters for building hybrid 
bonds. The wave function is automatically normalized within the QMC formalism. 

For the ratios V&new/^oid of the random walk and in the projection operator as well as in the parts contributing to 
the kinetic energy, the inverse matrices of the Slater matrices are conveniently determined numerically by the Sherman- 
Morrison- Woodbury formula. The periodic boundary conditions with respect to the wave functions are controlled by 
Green's relation - / vE^A* = / | V*| 2 . 
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The minimization over the huge parameter space is done as follows. Several hundreds of parameter sets in the region 
of the guessed minimum are chosen from a multidimensional Gaussian distribution. The total-energy expectation 
values weighted with their reciprocal statistical errors are fit into a multidimensional quadratic function. If the 
minimum of the fit lies within the region of the chosen statistical ensemble, a few parameter sets are selected along 
the minimum line of the fit with an increased number of simulation points. This calculation along the minimum line 
is statistically independent of the prior calculation with the large number of parameter sets and therefore not biased 
as it would be, if the absolute minima of several calculations were explicitly taken. 



III. RESULTS 



At first, a minimization for the ground-state energies of the gallium and the arsenic atom has been carried through. 
Table § gives an overview of the obtained quantities. The total energies are lower than the experimental values, 
which are the sums of the ionization potentials of the valence electrons. Only the valence electrons have been 
simulated in this approach. The differences are supposed to indicate the error from the frozen-core approximation 
using pseudopotentials, also apparent in the first ionization potentials. As expected, the £ parameters all result in 
1.0 from the minimization of the Hartree-Fock approach. Introducing the Jastrow factor we observe compression of 
the atomic orbitals even in the free atom. To deduce the constituents of the total energy, we made a Hartree-Fock 
simulation. The first ionization potentials are obtained from the approximation that the wave function is represented 
by the same parametrization scheme, i.e., by that of the neutral atom. 



TABLE I. Ground-state energies in eV of the gallium and arsenic atom from a variational QMC calculation: the first 
two lines are from Hartree-Fock calculations, the Hartree part of the electron-electron interaction is externally calculated, the 
exchange energy is the difference between the Hartree-Fock interaction energy and that Hartree part, the correlation energy is 
the difference between the total energy and the Hartree-Fock minimum total energy; the experimental total energy is the sum 
over the ionization potentials of the simulated electrons; the last line gives the first ionization potential. 



Energy 



Ga expt. 



Ga VQMC 



As expt. 



As VQMC 



Hartree-Fock, kinetic 
Hartree-Fock, potential 
Hartree, electron-electron 
Exchange, electron-electron 
Correlation, electron-electron 
Total 
Ionization 



-57.21 
5.999 



20.28 ±0.01 
-77.04 ±0.01 

42.08 ±0.04 
-16.21 ±0.04 
-1.250 ±0.008 
-58.000 ±0.003 
6.10 ±0.02 



-169.554 
9.81 



54.79 ±0.03 
-222.92 ±0.04 
134.3 ±0.7 
-37.3 ±0.8 
-2.00 ±0.02 
-170.134 ±0.007 
10.38 ±0.02 



FABLE II. Ground-state energies in eV per unit cell of gallium arsenide from a variational QMC calculation; i?km and Spot 
are produced by Hartree-Fock calculations, i? CO rr is the difference between E tot and the Hartree-Fock minimum total energy, 
Ecoh is the difference between £tot and the atomic energies, the experimental cohesive energy has been taken from Ref. [l7| the 
experimental total energy is the sum of the cohesive energy and the atomic total energies; do denotes the lattice constant of 
the zinc-blende crystal and B is the bulk modulus at temperature T = 0; LDA as well as the Becke and the Perdew gradient 
corrections (BP) are from density functional calculations of Ref. O. 



GaAs Expt. VQMC LDA BP 

~E~ 84.6 ±0.2 

E pot -311.3 ±0.3 

-Ecorr -6.42 ±0.2 

Etot -233.43 ±0.09 -233.04 ±0.08 

E coh -6.67 ±0.09 -4.9 ±0.2 -8.16 -6.45 

aoina.u. 10.6830 10.69 ±0.1 10.41 10.70 

Binkbar 756 786 ±100 
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Figure [l] shows the ground-state total energy of the gallium arsenide solid with 256 electrons as a function of 
the lattice constant, minimized with respect to all wave-function parameters. The curve is the quadratic function 
obtained from a weighted Gaussian fit of the data points. Table |0] shows the obtained quantities in comparison with 
the experimental ones and density functional calculations.EU The best agreement with the experiment is found for 
the total energy. The pseudopotentials give excellent results in this approach of parametrized linear combination of 
atomic orbitals wave functions even with the semilocal approximation. Differences in the cohesive energy originate 
from underestimating the atomic total energies of As and Ga owing to deficiencies of the pseudopotential in the atomic 
calculation. The lattice constant agrees very well with the experimental value, although this quantity statistically 
arises with a somewhat larger uncertainty than the total energy. The respective error given in Table [n] is an upper 
bound estimated by just considering the fit in Fig. |l]. 




The bulk modulus is obtained slightly too high, which indicates that lattice constants in the far neighborhood of 
the minimum are treated worse in the minimization procedure than those in the close neighborhood. This may be 
understood by observing that the minimum line is the result of the multidimensional quadratic fit, i.e., a straight line 
through the parameter space. In the far neighborhood of the global minimum this restriction could yield nonminimum 
parameters. The parts the total energy is composed of, as given in Table ||, are calculated in the minimum of the 
Hartree-Fock approach. The minimization was more difficult in this case, because the Coulomb cusp leads to a larger 
variance in the total energy. The difference between the results of the Hartree-Fock ansatz and the Jastrow-Slater 
ansatz, which defines the correlation energy, shows no statistically significant behavior in the simulated range of the 
lattice constant. Therefore the value of —6.42 eV per unit cell does not vary much near the ground-state density. 
The statistical data set for the correlation yields a straight regression line, the energy increasing with decreasing 
electron density with a slope of 0.9 ± 0.9 eV per r s . A measure for the ionicity is obtained by counting the electron 
visits to specific Voronoi polyhedra during a run. The polyhedra are constructed according to a partition of the 
crystal in which lattice sites are added at the empty tetrahedral positions and at the cube center, the so-called empty 
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spheres' positions in zinc-blende muffin-tin calculations. As a result the occupied As, occupied Ga, unoccupied As, 
and unoccupied Ga sites show charges of 3.9, 2.2, 0.9, and 1.0 electrons, respectively. Thus, both ions are positively 
charged the separated electrons extending into the interstitial space. The covalent bond charge was estimated by 
placing the centers of Voronoi polyhedra upon the bonds of the above lattice so that the polyhedra cover the bond 
region between sites that are occupied or empty. A slight charge accumulation appears on the occupied Ga-As bond 
of 0.13 electron per bond in addition to the charge which arises by the spherically symmetric electron distribution of 
the neutral atoms. 

Summarizing the obtained results, the good agreement with experiment puts the variational quantum Monte Carlo 
method into competition with current density functional techniques as far as the absolute ground-state quantities, 
calculated here, are concerned. The minimum principle guarantees an additional reliability absent in DFT, which 
in contrast bears an uncertainty because of the, in some respects, arbitrary choice of the density functional. As the 
latter is able to treat more complicated solids, the comparison with QMC calculations on the basis of simple systems 
may lead to a selection and adjustment of suitable functionals. 
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